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TECHNICAL NOTE 2820 


AN ANALYSIS OF THE ERRORS IN CURVE-FITTING PROBLEMS 
WITH AN APPLICATION TO THE CALCULATION OF 
STABILITY PARAMETERS FROM FLIGHT DATA 
By Marvin Shinbrot 

SUMMARY 


The problem of assessing the errors in the parameters obtained from 
a curve-fitting process is considered, and a scheme which may be applied 
toward the solution of such problems is obtained. This method is then 
specialized to the problem of finding the errors in the calculated 
stability parameters of an airplane, and an example is given. 


INTRODUCTION 


Curve-fitting procedures have found places in nearly all branches 
of engineering; in particular, the aeronautical engineer may apply these 
J methods to the calculation, from flight data, of the stability parameters 

of an airplane (references 1 and 2). Whether least squares or any of 
the profusion of graphical methods which exists is used for this curve- 
fitting process, questions of the errors in the calculated parameters 
are bound to arise. 

Although there is a considerable amount of literature on the subject 
of least squares and curve fitting, comparatively little is to be fo un d 
on the related subject of errors. What literature does exist (e.g., 
reference 3) attacks the problem from the point of view of statistics, 
arriving, finally, at a quantity called the variance. This quantity, 
while giving a satisfactory reply to the error question when applied to 
fitting a set of data to a straight line when only one measurement is 
subject to error, is far from adequate for other curve-fitting problems. 
One does not have to look far to find the reason for this; It is that no 
method has as yet been devised for calculating the variance when either 
the fitted curve is not linear or when more than one measured quantity 
is subject to error. This latter objection is not pertinent, perhaps, 
when the problem of calculation of stability parameters is considered, 

* for, although both the input (control-surface deflection) and the response 
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may be subject to error, it is frequently assumed that only the output 
is fallible. This approximation is particularly good when only free- 
oscillation data are analyzed. The first objection is, however, more 
serious, for, assuming even that all quantities remain within the so- 
called linear range, so that the response satisfies a linear differential 
equation, the response will certainly not be a linear function of the 
parameters. 

It is shown in the body of this report that the solution to the 
error problem depends upon a more or less arbitrary definition, the 
only criterion as to the choice of the definition being the usefulness 
of the solution to- which it leads. In this respect, it may be said of 
the method derived herein that it is not difficult to apply and appears 
to lead to reasonable values of the errors. A relationship between the 
error formula given herein and the classical formula for the variance 
(in the linear case when such a formula exists) is established. 


ANALYSIS 

Relation of the Problem to Aerodynamics 


Suppose one has a set of data which represents a time history of, 
for example, pitching velocity of an airplane in response to an elevator 
deflection &(t). If q(t) represents the pitching velocity, then, 
under certain simplifying assumptions, it is shown in reference 1 that 
q(t) satisfies the differential equation 


*Xb ia + 

dt 2 dt 


kq 


C, — + C 0 5 
1 dt 


( 1 ) 


where the constants b, k, Cj. and C 0 are functions of the stability 
derivatives of the airplane. 


It may be verified by differentiation that the following function 
represents the general solution of equation (l): 




q(t) = 


r A c i^i +c o r 

x,t 

P 4- 

[a , c i x 2 +c q rV x 2 T ^TW 

A^-f ■ — / e o^t ) U.T 

u 0 

C T 

X 2 -X .1 J Q 




where A, and A 2 are constants depending on the initial conditions q(o), 
(dq/dt) t: T 0 , and 5(o), while Xi and X 2 are "the roots of the character- 
istic equation 


* 


X 2 + bX + k = 0 



NACA TN 2820 


3 


It is often desirable to be able to ascertain the "best" values of the 
constants b and k (or X-i and \ 2 ) > C Q and Ci corresponding to the given 
time histories of q(t) and S(t). This problem is the subject of refer- 
ences 1 and 2 . 

Suppose then, that by some means these constants have been evaluated. 
There will be a certain "error" in these values, however, due to the 
experimental error in the data, the simplifying assumptions mentioned 
above, and other causes. The question of the magnitudes of these errors, 
given a certain error in the data, is considered in the present report. 


General Discussion of the Problem 


In order to state adequately the problem of errors, it is first 
necessary to give a precise statement of the curve-fitting problem. To 
this end, consider a physical quantity q e (t) which 1 b measured at t = t 0 , 
tf, ..., tjj, where t 0 <ti <...<tjj. It is assumed that various theo- 
retical considerations would indicate that q e ( t) should be one of the 
functions of the set q(t 1 x 1 ,x 2 , . . .jXjjj)) that is, there should exist 
values of the parameters x x , x 2 , ..., x^ such that 

q e (ti) = q(t ± , 

X 1> X 2> Xffl) ( 3) 

for all i =0, 1, ..., N. However, because of certain unknown errors 
in q 6 , equation ( 3 ) is not exactly satisfied for all i. It Is desira- 
ble then to find those values of the parameters which "most nearly" 

cause equation ( 3 ) to be satisfied. One means of doing this is to define 
the "best" values of x 1# . . . , x m as those values of the parameters 
which make 

N 

M = ^ [q(t^,x 1 ,x £ , . . . ,x m ) - qe(t^) ] (it-) 

i=o . 

a minimum. The process of minimizing M is called curve fitting by 
least squares, and this general problem is considered in reference 2 . 

As for the error problem, a careful study of the extant literature 
on statistics will lead one to the conclusion that the "true" value of a 
quantity, upon which the intuitive definition of error rests, has never 
been meaningfully defined. Clearly, however, the values obtained for the 
errors in a particular problem will depend upon this definition. When 
the usual probabilistic statement of the error problem is chosen, the 
quantity called the variance arises, along with _ the objections raised in 
the introduction which are concomitant with It. 
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Various nonprobabilistic statements of the problem are conceivable, 
and the statement given below has been chosen 1 from among them as a useful 
one. These nonprobabilistic statements are based on the concept of the 
"sensitivity" of the function q_(t,x 1 , . . . ,x m ) with respect to the parame- 
ters. That is, the .question is asked: If one of the parameters'. is 

changed slightly, does q also change only slightly, or does q change 
by a large amount? If a small change in a parameter ends in a large 
change in q, q is said to be sensitive to changes in that parameter, 
and it appears clear that any meaningful theory of errors should lead to 
relatively small error b in such a parameter. On the other hand, if q 
changes by only a small amount when a parameter is changed, q is insen-, 
sitive with. respect- to that parameter, and the theory should result in a 
large error. A quantitative discussion follows. 


Statement and Solution of the Problem 


In place of the elusive term "experimental error," we shall intro- 
duce the concept of "residual," defined by the equation 

e i - l( ti,X!,x 2 , . . . ,x m ) - qg(ti) (5) 

so that equation (4) may be written in the equivalent form 


N 



i=o 


Suppose the curve-fitting problem has been solved; that is, suppose 

values Xj/ 0 ), x 2 (°), •••> x m '°) of the parameters x 1 , x 2 , ..., x m , 
respectively, have been found which minimize M. Let 


ei (°) = 


= <1 




>*1 


(o) „ (o) 


,X 2 


> x m ' 


(o) 


- T e (ti) 


so that the minimum value of M is 


N 

M (o) . I 

i=o 

One possible statement of the error problem, which, with certain 
modifications, will be used in this report, is the following. Choose 
any set of numbers Ax x , Ax 2 , ..., A%, and let 

x k ^ ^ = x k (°) + Ax k (k=l,2, . . . ,m) 



* 
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The question then is: How large may the quantities | Ax&[ he so that 




+ . v (o) Y (o) (o) 

l ij x i > x 2 > * ■ • 


<Ih ( o) I (6) 


for all i = 0 , 1 , . . . , N? This maximum value of jAx^J may then be 
defined as the allowable error in x^. This problem can be solved, but 
various objections can be raised as to its significance as a statement 
of the error problem. In particular, it is easy to see that this estimate 
of the errors may be far too optimistic, for if only m out of the total 
of (N+l) values of ei( 0 / are zero, then the only way in which inequal- 
ity (6)' may be satisfied 1 for all i is for x^ 1 ) = x^f 0 ) for all k. 

That is, since the number m of parameters is usually far less than the 
number, (N+l) , of data points, we can say that even though most of the 
residuals be arbitrarily large, the result is that the error in the 
parameters is zero I 


This objection may be overcome by relaxing the requirement (6), 
substituting in its place a mean-square inequality obtained by squaring 
both sides of (6) and summing over i. The resulting inequality will 
then be required to hold instead of the inequality (6) . Symbolically, 
it will be required that 


N 


i=o 



• • > x m 


(1) 





(o) 





Thus, condition (6) is made to hold only in a mean-square sense over the 
whole range of t, while it may not, be true for some particular values 
of ti . An allowable error in x^ is then defined as any value of Ax^ 
for which inequality (7) holds. 


This last could reasonably be used as the definition of the error 
in xjj;. In the interest of ease in calculation, however, a Taylor's 
series expansion, with only the first-order terms being retained, will 
be used to linearize the problem. This may be done since neither the 
residuals nor the errors Ax^ may be too large; if they were, it 

may be said that there is something wrong with the theory which predicted 
that equation (3) will be approximately true, or with the experiment 
leading to the data q e (ti), or that the experiment has not been properly 
designed. In order to shorten the formulas, the following notation will 
be used: 1 

1 For the sake of clarity, it should be mentioned that two different types 
of subscript are used in this report, attached to the same quantity q. 
Of the experimental data q € , we have already spoken; in addition, the 
quantities q^- are now defined by equations (8) . There need be no 
confusion, however, for the subscript e will always be used to denote 
the data; other subscripts will be used as defined. 



6 


MCA TN 2820 


<l(ti) = q 




Q-k('ki) = “ * (k-lj2, • . • ,m) 


Sxt, 


We shall also write Aq(ti) in place of the difference 


( 8 ) 


t^,x^ - q Ft^jXi j • • 


so that the inequality (7) defining the errors in the x k becomes 

N 

l 2 


Y Aq(ti)J < 

i=o 


(7) 


If the errors in the parameters are not too large, it is clear that 
the following equation is approximately true; 


m 


Aq(ti) = Y q k (ti)Ax k 


k=i 


Utilizing this linearization process, inequality (7) becomes 

2 

< m(°) (9) 


The modifications of the definition of errors have finally been completed. 
We now define the error in x k as the largest value of [Ax k | for which 
inequality ( 9 ) holds, regardless of the values . of the other Ax's. 

In order to be able to draw some pictures, it will, now be assumed 
that m =-2. The generalization to larger m. will be presented after- 
wards. Thus, it is assumed that the function q is dependent on two 
parameters only: 


N 


i=o 


m 


Y Y qk(tj - )AXk 


k=i 


q = q(t i ,x 1 ,x a ) 



NACA TIT 2820 


7 


For "brevity in the formulas, 
define 


111 


*12 


0-22 


we shall write §£ 



\ 

[^(ti)] 2 


q.i( t i)q. 2 ( t i) ) 


[q. 2 (ti )] 2 


for Axjj. and shall 


( 10 ) 


The inequality (9) defines a certain domain in the (| 1J( S 2 ) plane whose 
boundary is given by the following equation: 


Q.11^1 2 + 2qx 2 iii 2 + Q. 22 I 2 2 



( 11 ) 


The graph of equation (ll) is either an ellipse, a parabola, or an hyper- 
bola. It is clear that if the definition of error which has been chosen 
is significant, equation (ll) must represent an ellipse, for if it did 
not, infinite errors would be obtained. It is just as well that the 
proof that equation (ll) does indeed represent an ellipse be given here, 
for it is very simple. Consider the discriminant 


In'! 


11^22 "*12 


of equation (ll), which becomes, using equations (10), 

N IT r N 

QnOea " <li2 2= £ M^i )! 2 £ h 2 {t±)] z - ^ q^tih^ti) 


-12 


1=0 


1=0 


L i=o 


The curve in question is an ellipse if and only if the discriminant is 
positive. However, the inequality 


IT 


N 


r n 


^ [^(ti) 3 2 £ [q 2 (ti)3 - £ qi(ti)q 2 (ti) 

i=o i=o _ i=o 


> 0 
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is precisely the "well-known Schwarz inequality. As is also well known, 
the equality cannot occur unless qj_(ti) is proportional to q 2 (ti) for 
all i, and as this implies that the parameters Xi and x 2 are not 
independent, it can he assumed, without loss of generality, that q x (ti) 
and q 2 (t^) are not proportional. This completes the proof. 

If the ellipse which represents equation (ll) should "be drawn, a 
graph similar to the one shown in sketch (a) would he obtained. 



Since the re are no first-order terms in equation (ll) , the ellipse must 
he sy mme tric with respect to the origin. Interpreting this curve, it 
may he said that any point which lies inside the ellipse has coordi- 
nates which define an allowable error in x x and x 2 . However, in the 
definition of the error in the parameters, it was said that the error 
in 2 % is the largest value of [Ax^l for which inequality ( 9 ) holds, 
regardless of the values of the other Ax’s. For this reason, in order 
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to define the error in the parameters "by means of the diagram, we must 
enclose the ellipse of sketch (a) in a rectangle, as in sketch (h) . 


^2 



Any point lying within this rectangle is considered as defining an 
allowable error. If H x and E 2 denote the errors (i.e., the maximum 
allowable errors) in x x and x 2 , respectively, then E x and are 
obtained as the maximum values of g 1 and | £ lying in the rectangle. 

If . m = 2 as In the above example, the errors can always be found 
graphically. However, If m>2, an analytical method must be used. 

Such a method will now be found for the case where m = 2, and the 
generalization will then be shown. 

It is clear that the sides of the rectangle drawn as solid lines 
in sketch (b) are the tangents to the ellipse at the points where 




= 0 



Calculating these derivatives from equation (ll), we obtain 
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_ _ lig ^1 + 9-aa £g 
d ^g = 9 U 5i + <l ls i 2 

> 

d ^g ^ In ^1 + 0-12 ^g 

d ^i q.ig 6i + <i 22 s 2 

4 

Setting these derivatives equal to zero results in the equations 


lig ^i + !gg ^g“ ® 


4n li+ CLig § 2 = 0 


( 12 ) 


These equations are not to he solved simultaneously, for they define 
errors in different parameters. Instead, the first of equations (12) 
is to he solved simultaneously with equation (11) to find Hi. the 
maximum value of Similarly, the second of equations (12) is to he 

solved along with equation (ll) to yield the error, H 2 , in x 2 . Per- 
forming these operations, one obtains the following expressions for Hi 
and *£ 2 : 



In the general case when m>2, the following equation occurs in 
place of equation (ll): 


where 


m 

I VA - m(o) 

j,k=l 


N 

Ijk = ^aj(ti)q k (ti) 

i=o 


( 1*0 


( 15 ) 
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and the quantities CLjfai), qk(t£) are defined in equation (8). If Hh. 
denotes the maxi m u m , value of satisfying equation (l4) , so t ha t Hh. 

is the error in , the problem is to find Hu. It may be solved by 
setting the (m-l) derivatives of the form 


( « 3 = 1 , 2 , • • . ,h-l,h+l, • . . ,m) 

equal to zero and solving the resulting equations along with equation (lit-) 
fo*" £b* not go through this process each time the errors in 

a problem are to be calculated, for it can be done once and for al 1 as 
follows: Differentiating equation (l4), it may be seen that 



Setting these derivatives equal to zero, we obtain the following (m-l) 
equations which are to be solved simultaneously along with equation (1^) 
for 


Equation (1^) 


m 

y q^jj = 0, j=l,2, ... ,h-l,h+l, ... ,m 

sA 

can be written in the form 


(16) 



/^ia 


/ 





( 1*0 


Taking equation (l6) into account, it may be seen that the parenthesized 
quantity is zero except when j=h, and so we obtain the equation 


m 


Yj ^lk^k = 

k=x 


Jf(f) 

^h 


(IT) 


This equation and equation (16) are m linear equations to be solved 
for 1^. This may be done, in general, by means of Cramer’s rule (the 
method of determinants) . 
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Let 


Dh 




^•n ^ 12 ’ * ’ 9-im 




D = 

9 21 9-22 " ’ * 9pm 





9mi 9.JH2 * ' * 1mm 



be 

the minor of q^ in D, so that 


111 

9-12 

<1 i(h-i) ^i(h+i) 

• • • 

^im 

9-21 

q 22 ... 

q a(h-i) C1 2(h+i) 

• • • 

9-2m 


^(h-i) 2 * * * 

q -(h-i)(h-i) q -(h-l)(h+i) •** 

q -(h-i)m 

^(h+i) l 

• • • • 

q (h+l)2 

q (h+i) (h-i)*(h+i) (h+i) *•* 

^(h+i)m 

* • • 

^mi 

^2 

9m(h-i) 9 m (h+i) 

• • « 

1mm 


Then, "by virtue of equations (l6) and (17), 

m(°) % 


= 


(18) 


(19) 


ty 


D 


That is. 


~h 



Equation (20) is the desired formula for the error in 


( 20 ) 


Connection of formula (20) -with the variance concept .- Although the 
formula (20) for the errors was obtained without the aid of statistical 
notions, it may now he shown that a relationship exists between the 
parametric errors as defined herein and the concept of variance. 
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Suppose that the given approximating function q(t) were linear in 
the parameters, so that 

q(ti).« 1 1 (ti)x 1 + <i 2 (t i )x 2 + ... + ^(ti)^ 


where the functions qj (t-j_) are independent of the para m eters 
Differentiating this last equation, ve have that 


- = qk(ti) (k=l, ... , m) 


and, utilizing equation (15), that 


* ^ <lj(ti)q k (t i ) 
i=o 

In reference 3j the familiar normal equations obtained from the least- 
squares process are shown to he 



where 2 


m 

£ WJ = (21) 

J-f 

N 

q e ,k = ^ qe(ti)q k (ti) . 
i=o 


.and q e (tj_) represents the given data. Let D denote the determinant 
of the coefficients In equation (21), and let Eh he the minor of q^ 
In D, so that D and are given by equations (l8) and (19) , respec- 
tively. Then, if 0 2 denotes the mean-square error in the data, while 
b k 2 denotes the mean-square error to he expected in the parameter x k 
(i.e., b 2 and 0 k 2 are the variances of the data and of x^, respec- 
tively), it is shown in reference 3 that if q(t^) is a linear function 
of the parameters. 



2 The co mma is placed after the subscript e to indicate that it is a 
different sort of symbol from the other subscripts. Thus, j and k are 
to vary, assuming the values 1, . . . , m; however, the subscript e is 
used merely to indicate that the experimental data q e (tj_) is used in 
the definition of q e ^ k * 
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Noting now that a 2 affords an estimate for the error in the data when. 
the statistical conception of the error problem is considered, while 
M^°) occupies a similar position in the present theory, a relationship 
between the variance in a parameter and what has herein been called the 
allowable error in that parameter can be established. Comparing equa- 
tions (20) and (22), it -may be seen that the ratio of the variance in a 
parameter to the statistical measure of the error in the .data is the 
same as the ratio of the error in a parameter to the nonstatistical 
measure of the error in the data. In symbols, 



(23) 


This is a remarkable fact that in formulating a theory of errors which 
abstains from the use of probability theory, a notion of error has still 
been defined which bears as intimate a relation to the statistical ideas 
as that described by equation (23). Of course, this relationship has 
been proved only in the case of linear curve -fitting problems. From 
this, however, it might appear reasonable that even for general nonlinear 
problems equations (20) and (23) (and, therefore, (22)) afford an expres- 
sion for the variance. This leaves an interesting problem for future 
research. 

There is an important conclusion to be drawn from equation (23) • 

If, as before, N + 1 is the number of data points, then a 2 and m(°) 
are related (reference 3) hy the equation 



N+l 


Thus, from equation (23 ) > 


H t 2 = (N+l)a h 2 - 

which implies that is larger than 1 %. This result should not be 

too surprising, for It will be recalled that c^ 2 is a measure of the 
probable error in x^, While we have tried, ^in to define some sort 

of maximum error. — ■ 


APPLICATION TO AN EXAMPLE 


Referring to equation (2), if there Is some value to, of time such 
that 5(t) =0 for all t>t a (i.e., if a pulse elevator input has been 
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applied to the airplane), it may he seen that for t>t a , the pitching 
velocity is a sum of two exponentials 


ti 

1 e (ti) 

0.4 

0.224 

•5 

.120 

.6 

.020 

.7 

-.057 

.8 

-.112 

• 9 

-.148 

1.0 

-.160 

l.l 

-.150 

1.2 

-.127 

1.3 

-.097 

1.4 

-.062 

1.5 

-.032 

1.6 

-.005 

1.7 

.017 

1.8 

.030 

1-9 

.036 

2.0 

.035 

2.1 

.032 

2.2 

.027 

2.3 

.020 

2.4 

.015 

2.5 

.011 

2.6 

• 008 

2.7 

.005 

2.8 

.003 

2.9 

.001 

3.0 

0 

3.1 

0 

3.2 

0 


q = + Bse^^t 

where B 1 ,B 2 ,A 1 and X 2 are constants depending on the 
stability derivatives of the airplane. 

In general, such data are oscillatory; that 
is, 1 1 ,1 2 ,Bi and B 2 are complex. Let = l -f-z’i, 

(g+g'i), i 2 = - 1. Since q is real, \ z = l-V i, 
B 2 - | (0-0*i), and 


q = e^(0 cos Z’p - 0' sin Z’t) 


(24) 


The data given on this page are actual flight data 
which represent the pitching velocity of a test airplane 
in response to an elevator deflection which was zero for 
all t >0.4 (t * 0 is taken at the beginning of the pulse). 
These data are fitted to a sum of exponentials in refer- 
ence 2, where the results 


Z = -I.366 0 - 0.614 

z* = 3.071 p* = -0.208 , 

m(°) = 0.000895 


(25) 


were obtained. The errors in these parameters will now be 
calculated, utilizing equation (20) and substituting Z , Z*, 0, 
and g’ for x 1 ,x 2 ,x a , and x 4 , respectively. 


Calculation of the errors in l, V, g, and P 1 . - Taking derivatives 
from equation (24), one may see that 


= te^(g cos Z't - g’ sin Z’t) 


8q 

dZ 

— = -te^(g sin Z’t + g’ cos Z’t) 
81 ’ 


h-. e it 

8g 


cos 


Z't 


— - -eZt s in z r t 

8g» 




(26) 
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and 


Identifying 2 with, the parameter 
P* with x 4 , it may he seen that q (tj.) 

column 11 of table I and q 2 (ti) is given by minus column* 15, while 
q. 3 Jti; and q. 4 (ti) are given, respectively, by column 6 and by minus 
column 7. Therefore, referring to equations (15), and letting circled 
numbers refer to columns in table I, we have that 


with x 2 , p with x a , 

Sq(ti)l 4 

= — — — J is given in 


111 - Z @ 2 = 0.168 

'll 2 = <l 2 i = -Z @ x 

Q 

* 0.006 

<*13 = 'l 31 = £ O X 

@ 

* 0.228 

*14 " *41 = -2 © X 

© 

= -0.091 

I22 = 2 © 2 ■ 0.212 
*23 = I32 " - s © X 

© 

= 0.117 

*24 " *42 = 2 © X 

© 

= O.389 

*33 = z © 2 = 0.415 

*34 = *43 = -2 ® x 


> 0.051 

2 

*44 " Z © - 0.985 


Inserting these numbers into the expression (18) for D, we obtain 


0.168 

0.006 

0.228 

-0.091 

.006 

.212 

.117 

.389 

.228 

.117 

.415 

.051 

-.091 

.389 

.051 

.985 


Furthermore, 


0.000343 


V-’ 

0.212 

.117 

0.117 

.415 

. 0.389 
■ .051 

= 0.0145 


.389 

.051 

’ .985 

d 2 = 

0.168 

.228 

0.228 

.415 

-O.091 

.051 

= 0.0115 


-.091 

.051 

.985 
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0.168 

0.006 

-0.091 


d 3 = 

.006 

.212 

.389 

= 0.00744 

-.091 

.389 

.985 



0.168 

0.006 

0.228 


d 4 = 

.006 

.212 

.117 

* 0.00177 


.228 

.117 

.415 



Substituting these numbers and the value of m(°) given in equations (25) 
into equation (20), we obtain 


|AZ [= 0.194 
j Az* | = .173 „ 
|AP (= -139 
j AP * j = .068 


(27) 


where AZ denotes the error in z, AZ' in Z', etc. On a percentage 
basis, this implies that Z, Z 1 , p, and p ' , respectively, are known to 
within about 14, 6, 23, and 33 percent. 


Calculation of the error in the stability parameters .- It has been 
shown how the errors in the parameters i' t l x , p, and P* may be found. 
However, the problem of interest to the aerodynamicist is the calculation 
of the errors in the stability parameters b, k, C^, and C Q of equa- 
tion (l). 'While q(t) can be written directly as a function of b, k, Ci, 
and C 0 instead of the parameters Z, Z*, P, and p T and the method 
described above applied. It is believed that the following general con- 
siderations simplify the calculations. 

Suppose, as before, that q(t) is a function of the m param- 
eters x k , k=l, ...,m, and that the errors Ax^ in the xk have been 
found. Suppose further that there are n other parameters y 1 , . . . ,y n , 
each of which is a function of xi, x 2 , ..., and xja* Thus, 


y^ = yj (x^, • * • , x m) , j-l#...^a 


The error Ayj in yj may then be estimated from the formula 


hi- J 

k=i 


Syj 

5x k 


Ax k 


(28) 


It should be noted that equation (28) gives a quite pessimistic value 
of the errors. 
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From the definition of l and 2', we may now write 


and so, 


Also, 


b 

k 


-22 

l 2 + 2* 2 


}■ 


8b 2 8b m 8b ^ 8b 

bl ~ ’ 82 ' 8p 8p* 


8k 

8k 

— a 22 
82 

8p 

^ = 22' 

8k 

82 ' 

8p» 


■ 0 


= 0 


(29) 


Thus, applying equation (28) using the error values (27), we see that 

|Ab | = 2 | A 2 |= 0.388 

|Afc | = 2 | 2A2 |+2 | I'AV | = 1.59 

and since, from equations ( 29 ), 

b = 2.732 | 
k = 11.30 J 

it follows that b and k are known to within about l4 percent. 

This same scheme will now be applied to the computation of the 
errors in C x and C c . It will be noticed that since Ci and 0 o are 
quite complicated functions of 2 and 2', our -task will be considerably 
more difficult than it was when computing |Ab| and |Ak|. 

As before, suppose the input 5(t) is zero for all t greater 
than ta. Make the following definitions: 
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/ la 

e"^ T 5(T) cos Z'TdT 
*i 

/ ^a 

e“^ T 5 (t) sin Z'TdT 

H 

ft* 

Rj_ = / Te"^ T 5(T) cos Z'TdT 

ft* 

Ri' = / Te"^ T 6(T) sin Z'TdT 


so that 


dSi dSi' 


-Ri 


SSi dsi» 


— Rj ' 


Also define 


Oqj = e^CSi* cos Z'ti - Si sin Z'li) 

a li = -e^CCZ'Si-ZSi') cos Z'ti + (zSi+Z'Si’) sin Z’ti] 


a i ~ l 1 [ 0. 6 (li ) " e cos " P* sin Z'li) ) 


a oo = 7 


-l 


void'll ao = > ccoicci 


■I 


a io - 7 a n 


it 

■I 


it 

-I 


a ii a i 


With these definitions f it is shown in reference 2 that C 0 and C x 
satisfy the equations 


aoo C 0 + a Q1 Cj_ — &q 


a lo Cq + a i; i_ Cj_ = aj 
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Therefore, if p denotes any one of the four parameters l, l 1 , 3 , 3 r , 
the derivatives 8c 0 /8p and 8c x /8p can he found from the equations 


a oo 


l io 


dC 0 ^ 

8 C 1 

8 a 0 

_ ^ a OO 

8 aoi 

-r— + Spi 

op 

8 p 

Sp 

0 dp . 

1 Sp 

S0 o 

8 Ci 

8 a! 

„ Sa io 

8 a n 

IT + 811 

8 p 

8 p 

c ° a P 

1 

0 

H 

& 


( 35 ) 


The problem before us is the calculation of these derivatives. 

From the definitions (32) and from equations (31), it follows that 

= e Zti [(tiSi« - Ri>) cos Z’ti - (tiSi - Ri) sin Z'ti] 
o Z 

3 a ii 

-g— = tiaji + otQi + e 1 [( Z'Ri - ZRt') cos Z'ti + (zR t + z'Ri') sin Z'ti]} 

8a j_ 


8 Z 


= - Z , t i e iti (p cos Z’ti - 3» sin z’ti) 


~~ - eUl [(Ri “ ttSi) cos Z’ti + (Ri» - t^*) sin z’tq] 


8z 

= -e Ui {[Si + j(tiSi - Ri) + Z»(tiSi» - Ri*)I cos Z't ± + 
[Si 1 - ^'('tiSi - Rq) + z(tiSi f - Ri 1 )] s i n Z'ti 


8aj OU 7-fcj , 

— = — + Z’tie 1 (p sin Z’ti + 3 ' cos Z'ti) 


> 


8a, 


oi 


8p 

8aii 

83 

8 ai 

83 

dapj 

83' 

8aii 

83’ 

8ai 

83' 


= 0 


= 0 


= -Z'e^ i cos Z'ti 


= 0 


Z'e^ sin Z'ti 
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The quantities e" lt 8(t) cos l't f e" lt B(t) sin I’t, te“ Zt 8(t) cosl’t, 
te“^8(t) sin Z’t must he integrated to find Si, Si’, Rj_,andRi’. 

These quantities are calculated in table II and are plotted versus time 
in figure 1. They were integrated by means of a planimeter to obtain 
the values displayed in table II. Referring to table II and equa- 
tions (33), it may be seen that 


a-oo » 2 @ 2 = 0.000399 

®oi “ a 10 - S @ _ X @ 
2 

a ll=2© = 0.02feL 


-0.00292 


a o =2 @ x @ = -0.0615 
a i = 2 © X © = 0 : k 9 k 


and that 


8a 


00 


= 22 


8z 

8a 01 8a 10 


8l 

^ a n 

81 

8a 0 

8i 

8a x 

ITT 


8z 


2 2 


= 2 


= 2 


= 2 


X 


= -0.000152 


X 


+ 2(3(3 x 


= 0.001388 


= -0.0l26i+ 


+ 2(360 x 


+ 2Rd) x 


0.00772 


= -0.0851 


—22. =22(0) x © = 0.000232 
8 1 1 

—I = »2 @x©+ 2 @x©= -0.001571 

al l oV 

= 2 2 @ X © = 0.01226 
^ = 2@ X @ +2 @ X © = -0.0381 

= 2 © x © + 2 © x © = 0.280 
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8a 00 

= 0 



Sa 00 

= 0 


- 

c> a oi 

9a io 

K ■ ■ ' mm 

ss 0 


8a 01 

8a 10 

— n 



8(3 



83* 

83' 

A 



Sa n 

- 0 



^ a ll 

— n 






83' 

— u 



8a 0 

-s @ 

x © 

■ -0.10*4- 

8a 0 

83* 

- z @ 

X © 

= 0.0163 

8a x 

= s© 

X © 

= 0.899 

8a j. 
83' 

= S® 

X © 

= -0.22*4- 


Substituting these numbers in the appropriate places in eq.ua- 
tions (3*0 and (35), it follows that 

C 0 = -37.90 
C x = 15.88 

while 


^2 . 6.65 

8 1 

8c x 

Hi m 7 - 79 

8c 0 

r- 2 - = -27.38 
8 1 ' 

8c, 

= -2.25 
8z» 

8c 0 

~ - 108.93 
83 

8c, 

~ = 50.50 

83 

8c 0 

r- 2 « -239.79 
83* 

8c, 

r-i = -38.35 
83* 


Using equations (28) and (27), we obtain, finally, 

}ACo!= 37.5 
\ac 1 \= 11.5 

Therefore, C3. is known to within about 72 percent; C 0 is only known to 
within 99 percent. 



MCA TN 2820 


23 


The large errors in some of the parameters are worthy of 
special comment. It should be noted first of all that the particularly 
large error in C 0 is completely consistent with what has been found 
for some time by empirical means. It has been found that upon repeat- 
ing a set of flight records, making the two sets identical as far as 
that is possible, two entirely different values of C Q have often been 
obtained. 

Secondly, it will be recalled that in the derivation of the formula 
for the errors it was required that the errors be small. j this is cer- 
tainly not the case for some of the errors. For this reason, the error 
of 37-5 in C 0 particularly may not be considered as definitive. The 
theory developed herein fails for such a large error, and the actual 
error in Cq may be much larger than that calculated. However, this is 
not serious, for an error as large as 99 percent renders the calculated 
parameter value meaningless in any case, and it matters not at all 
whether the error is 99 or 199 percent. 

Finally and most important is. the following conclusion regarding 
the entire philosophy of calculation of stability parameters from pitch- 
ing velocity data alone in response to an elevator pulse: This experi- 

ment is ill designed for the calculation 3 of Co. In. general, when a 
pulse is applied to produce a set of pitching velocity data, no more 
should really be expected from the analysis than the period and the 
damping parameters k and b. If, in addition, the analysis results in 
a value of Ci whose error is within reasonable limits, this should be 
regarded as fortuitous. It should be stated that no example has yet 
been found for which C D nay he calculated at all accurately. 

Such a negative comment as the preceding deserves a remark on the 
possibility of the most accurate calculation of Co* If one were to 
measure a step response rather than a pulse response, the analysis of 
these data would surely lead to more reliable values of Co* The 
reason for this lies in the easily proved proportionality of the steady- 
state value of q.(t) with C 0 . 


3 This conclusion is drawn for the parameter C Q alone and not also 
for Ci, as would seem to be indicated by the error of J2 percent 
in Ci, primarily because the error in Ci is not always as large as 
this. With some data which have been analyzed, the error in C x has ' 
been far less, of the order of 30 percent. 
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CONCLUDING REMARKS 


A formula has "been given which may he used to find the errors in 
the parameters obtained from a curve -fitting process. The method of 
derivation did not use the concepts of probability theory, since the 
latter ideas lead to quantities for which there is no known method of 
solution. However, the result that the formula obtained bears a close 
relationship to the classical probabilistic formula in case the curve- 
fitting problem is linear (the only case for which such a formula has _ 
been derived and proved valid) is proved. 

As may be seen by studying the example given, the ease of appli- 
cation of the method is directly dependent on the simplicity of the 
form of the function which has been fitted. Thus, if b, k, Cx, and C 0 
are the aerodynamic parameters occurring in. the differential equation (l) 
the errors in b and k may be calculated fairly rapidly; the errors 
in Cx and C G require more time. However, if the method of least 
squares (reference 2) is used for the curve -fitting problem, a great 
many of the computations which are needed to find the errors will have 
already been performed in the process of finding the parameters. 

The method actually Weights favorably the method of least squares 
in another mann er. The errors, as given by equation (20), are propor- 
tional to the square root of the sum of the squares of the residuals, 
which quantity is minimized by the least-squares process. However, it 
is not believed that this weighting is a serious limitation, for no 
matter what means are used for fitting the curve, if this latter method 
is to have significance and is to lead to a good fit of the data, the 
sum of the squares of the residuals will also be near its minimum. 

Finally, it has been shown that pitching velocity data in response 
to a pulse are not alone adequate to compute the parameter C 0 (and, to 
a lesser extent, Cx) occurring in the differential equation of motion. 
Usually, all that may be obtained from such data are the parameters b 
and k, which determine the damping and the period of the oscillation. 

In some cases, C x nay also be obtained with reasonable accuracy. 

Ames Aeronautical Laboratory 

National Advisory Committee for Aeronautics 
Moffett Field, Calif., July 30 , 1952 
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TABLE I.- THE ERRORS IE l, l' , p MD p* 


1 

2 

3 

4 

5 

l«J 

I 

8 

9 

10 

n 

12 

13 

14 

15 

Raw 

■^i 


cos Z't^ 

sin Z' tjL 

( Dx © © x © 

©X 0 

©X 0’ 

®-(B 

©X© 

©X 0 

©X 0* 


1 

0.4 

0.579 

0.336 

0.942 

0.194 

0.545 

0.119 

-0.U4 

0.233 

0.093 

0.335 

-0.041 

0.294 

0.118 

2 

.5 

.505 

.035 

.999 

.018 

• 505 

.011 

-.105 

.116 

.058 

• 310 

-.004 

.306 

.153 

3 

.6 

.441 

-.269 

.963 

-. 118 

.424 

-.073 

-.088 

.016 

.009 

.261 

.025 

.285 

.171 

4 

• 7 

.384 

-.547 

.837 

-.210 

.322 

-.129 

-.067 

-.062 

-.044 

.198 

.044 

.241 

.169 

5 

.8 

.335 

-.775 

.632 

-.260 

.212 

-.160 

-.044 

-.115 

-.092 

.130 

.054 

.184 

.147 

6 

• 9 

.292 

-.930 

.369 

-.272 

.108 

-.167 

-.022 

-.144 

-.130 

.066 

.057 

.123 

.111 

7 

1.0 

.255 

-.998 

.070 

-.254 

.018 

-.156 

-.004 

-.152 

-.152 

.011 

.053 

.064 

.064 

8 

l.l 

.223 

-.972 

-.235 

-.216 

-.052 

-.133 

.011 

-.144' 

-.158 

-.032 

.045 

.013 

.014 

9 

1.2 

.194 

-.856 

-.518 

-.166 

-.100 

-.102 

.021 

-.123 

-.148 

-.062 

.035 

-.027 

-.033 

10 

1.3 

.169 

-.659 

-.752 

-.112 

-.127 

-.064 

-027 

-.095 

-.124 

-.078 

.023 

-.055 

-.071 

11 

1.4 

.148 

-.401 

-.916 

-.059 

-.135 

-.036 

.028 

-.065 

-.090 

-.083 

.012 

-.071 

-.099 

12 

1.5 

.129 

-.105 

-.995 

-.014 

-.128 

-.008 

.027 

-.035 

-053 

-.079 

.003 

-.076 

-.n 4 

13 

1.6 

.112 

.200 

-.980 

.023 

-.110 

.014 

.023 

-.009 

-.015 

-.068 

-.005 

-.072 

-.116 

l4 

1.7 

.098 

.487 

-.873 

.048 

-.086 

.029 

.018 

.oil 

.020 

-.053 

-.010 

-.Q63 

-.106 

15 

1.8 

.086 

.728 

-.685 

.062 

-.059 

.038 

.012 

.026 

.047 

-.036 

-.013 

-.049 

-.088 

16 

1.9 

.075 

.902 

-.433 

.067 

-.032 

.o 4 l 

.007 

.035 

.066 

-.020 

-. 6 l 4 

-.034 

-.064 

17 

2.0 

.065 

• 990 

-.l 4 o 

.064 

-.009 

.o 4 o 

.002 

.038 

.076 

-.006 

-.013 

-.019 

-.038 

18 

2.1 

.057 

.987 

.166 

.056 

.009 

.034 

-.002 

.036 

.076 

.006 

-.012 

-.006 

-.012 

19 

2.2 

.050 

.890 

.456 

.044 

.023 

.027 

-.005 

.032 

.070 

.014, 

-.Q09, 

.005 

.qio 

20 

2.3 

•043 

.710 

.704 

.031 

.030 

.CCL9 

-.CI06 

.025 

.058 

.019' 

-.006 

.012 

.028 

21 

2.4 

.038 

.464 

.886 

.018 

.033 

.011 

-.007 

.018 

.043 

.021 

-.oo4 

.017 

• 040 

22 

2.5 

.033 

.175 

.985 

.006 

.032 

.004 

-.007 

.010 

.026 

.020 

-.001 

.019 

.047 

23 

2.6 

.029 

-.131 

.991 

-.004 

.028 

-.002 

-.006 

.oo 4 

.009 

.017 

.001 

.018 

.047 

24 

2.7 

.025 

-.425 

.905 

-.011 

.023 

-.007 

-.005 

-.002 

-.005 

.014 

.002 

.016 

.044 

25 

2.8 

.022 

-.678 

.735 

-.015 

.016 

-.009 

-.003 

-.006 

-.016 

.010 

.003 

.013 

.036 

26 

2.9 

.019 

-.869 

•495 

-.017 

.009 

-.010 

-.002 

-.008 

-.024 

.006 

.003 

.009 

.027 

27 

3-0 

.017 

-.978 

.209 

-.016 

.004 

-.010 

-.001 

-.009 

-.027 

.002 

.003 

.006 

.017 

28 


.015 

-.995 

-.096 

-.014 

-.001 

-.009 

0 ' 

-.009 

-.028 

-.001 

.003 

.002 

.007 

29 


.013 

-.920 

-•393 

-.012 

-.005 

-.007 

.001 

-.008 

-.026 

-.003 

.002 

-.001 

-.002 


= 0.168 

- -0.0p6 
= 0.228 

- 0.091 
= 0.212 


£ ©x © 

*©* © 
£( 6 ) 

S®x© 

*©* 


—0.117 
= 0.389 
= 0.415 
- -0.051 
= 0.985 
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TABLE II.- THE ERRORS JN Ci ARD C Q 


rn 

m 

m 

m 


6 

T 

a 

9 

m 

11 

12 

13 

14 

15 

16 

n 

ra 

19 

20 

21 




Ig 


•In Z*ti 





Si 

Si* 

®x® 


Ri 

Hi* 

m 

a 


@x@ 

J*ol - _ 

© - @ 

K 



llwl 

1.000 

0 

1.000 

m 

ran 

RH 

0.02172 

0.01490 

0 

0 

0.004186 

0.003234 

1.000 

0 

0.01498 

0 

0.01498 

K 

.05 

.005 

1.071 

.968 

.153 

1.058 

.16* 

.005 

.001 

.02167 

.01497 

.00025 

.00005 

.004184 

.003934 

.9e3 

.143 

.01382 

.00310 

.01072 

■ 

.10 

.0*5 

1.1M 

.953 

.302 

1.093 

.3^7 

.049 

.016 

.02CA 

.01468 

.00490 

.00160 

.004074 

.003210 

.032 

.264 


.00540 

.00681 


.1? 

.105 

1.227 

.896 

.4*5 

1.09? 

•546 

.115 

.057 

.01636 

.01278 

.01725 

►00055 

.003X4 

.002978 

.730 

.363 

.00933 

.00595 

.00338 

m 

.20 

.125 

1.3U 

.817 

.576 

1.07V 

.757 

.13V 

.095 

.00996 

.009CC 

.02680 

.01900 

.002410 

.002298 

.622 

.438 

.00560 

.00435 

.00124 

m 

.25 

.100 

1.407 

.719 

.695 

1.012 

• 977 

.101 

.098 

.003 60 

.00410 

.02525 

.02450 

-001050 

.oou4a 

.511 

.494 

.00210 

.00178 

.00032 

■ 

.* 

.030 

1.50T 

.605 

.796 

■9U 

1.200 

.027 

.036 

.00040 

.00076 

.00810 

.01080 

.000150 

.000126 

.401 

.526 

.00031 

.00021 

.00010 

E 

• 3? 

0 

1.613 


.660 

.768 

1.419 

0 

0 

0 

0 

0 

0 

0 

0 

.295 

.546 

0 

0 

0 

E 

.kc 

0 

1.727 

.336 

.942 

.560 

1.6e7 

« 

0 

0 

0 

0 

0 

0 

0 

.195 

,545 

0 

0 

0 


r 

22 

23 

24 

25 

26 

27 

26 

29 

*> 

MU 

32 

33 

■a 

33 

Ell 

EH 

38 1 


@x r 

@x l 

©-© 

@X @ 

@X t 

@xi> 


©X@ 

*11 - 
-[© + © 

EH 

©Xf 




°i - 

©X@ 


1 

E 

0.06670 

.06655 

.06277 

.05030 
•0*59 
.01106 
.00123 
0 - 
0 

0.02046 

-.02045 

-.03305 

-.01746 

-.01229 

-.00560 

-.00104 

0 

0 

O.ofrP-6 

.08700 

.08082 

.06776 

.04288 

.01666 

.00227 

0 

0 

0.06720 

.03021 

.06897 

.04949 

.02688 

.00853 

.00088 

0 

0 

,0.02967 

-.02960 

-.00792 

-.00238 

-.01361 

-.00492 

-.00055 

0 

0.C4500 

.04597 

.04508 

.03925 

.02764 

.01259 

.00233 

0 

0 

0.01633 

.01637 

.01716 

.01687 

.01403 

.30767 

.00178 

0 

0 

0 

.00234 

.00453 

.00612 

.00615 

.00379 

.00094 

0 

0 

-0.06720 

-.08255 

-.07350 

-.05561 

-.03263 

-.01232 

-.00182 

0 

0 

0.6l4 
.567 
.511 
.448 
.362 
• 314 
.246 
.181 
.120 

0 

-.0* 

-.055 

-.076 

-.091 

-.103 

-.110 

-.114 

-.113 

o.6l4 

.597 

.566 

.524 

.473 

.417 

.356 

.295 

.233 

0 

.006 

.095 

.200 

.310 

.360 

.320 

.272 

.224 

-0.614 

-589 

-.471 

-.324 

-.163 

-.057 

-.036 

-.023 

-.009 

-1*806 

-1.809 

-1.446 

-.995 

-.501 

-.175 

-.111 

-.071 

-.028 

.0007*? 

.001466 

.001917 

.001800 

.001025 

.000228 

E 

-0.003234 

-.002485 

-.001742 

-.001061 

-.000496 

-.000123 

.000100 

0 

0 

■ 

39 

40 

41 

42 

43 


45 

46 

47 

48 

49 

50 

51 

52 

■9 



8 

@X@ 

©*@ 

©-© 

©x© 

*01 _ 
31 7 
©-© 



©X !• 

©x 1 

m 

0X0 


@xl> 

©+@ 

©X© 

n 

m 

1 

2 

3 

4 

5 

6 

7 

8 
9 

0.003234 

-.002294 

-.001449 

-.000775 

-.000310 

-.000063 

.00004Q 

0 

0 

0 

.001064 

.002CU 

.002457 

•00199C 

.000900 

.000120 

0 

0 

■0.004186 

-.003100 

-.0020* 

-.001047 

-.000418 

-.000150 

-.0000* 

0 

0 

0 

-.000443 

-.000536 

-.000380 

-.000183 

-.000074 

-.000016 

0 

0 

-0.003234 

-.001851 

-.000913 

-.000395 
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